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Abstract 



We study the perturbative QCD series for the hadronic width of the Z 
boson. We sum a class of large "vr^ terms" and reorganize the series so as to 
minimize "renormalon" effects. We also consider the renormalization scheme- 
scale ambiguity of the perturbative results. We find that, with three nontrivial 
known terms in the perturbative expansion, the treatment of the vr^ terms is 
quite important, while renormalon effects are less important. The measured 
hadronic width of the Z is often used to determine the value of as{M'^). A 
standard method is to use the perturbative expansion for the width truncated 
at order in the MS scheme with scale fi = Mz- We estimate that the 
determined value of as{M^) should be increased by 0.6% compared to the 
value extracted with this standard method. After this adjustment for vr^ and 
renormalon effects, we estimate that the uncertainty in as{M^) arising from 
QCD theory is about 0.4%. This is, of course, much less than the experimental 
uncertainty of about 5%. 
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I. INTRODUCTION 



The width for Z hadrons is conventionally described by the ratio R of this width to 
the width for Z e^e~ . The Z boson need not be on-shell: for theoretical purposes, we can 
consider i? as a function -R(s) of the cm. energy s of the e^e~ annihilation that produces 
the Z. Then the measured R is R{M'^). One way of measuring the strong coupling as is to 
compare theory and experiment for R^M"^). The purpose of this paper is to discuss some 
aspects of the theoretical evaluation of R{s): the effect of "vr^ terms" and "renormalons" on 
the determination of R from the calculated terms in its perturbative expansion in powers of 
as- Our goal is to suggest ways of evaluating R{s) as precisely as possible from the knowledge 
of the first three terms in its perturbative expansion and then to estimate the theoretical 
error in this evaluation. We pose the question of whether as(M|) could be extracted at a 
precision of a few parts per mill from i?(M|) in the hypothetical case that infinitely accurate 
data were available and uncertainties in the electroweak part of the calculation were zero. 

We will conclude that a QCD theoretical error on q;s(M|) of about four parts per mill is 
possible if one understands this as a one a error estimate: the QCD error is probably about 
this size. An estimate of the QCD theoretical error at the 95% confidence level would be 
quite a lot larger because it should include the possibility that certain hypotheses - guesses 
really - about the behavior of the perturbative expansion are simply wrong. We will try 
to make clear the nature of the required hypotheses and let the reader form his or her own 
judgment. 

In this paper, we adopt a simplified theoretical framework so that we can concentrate on 
the QCD effects. We consider Z — > hadrons at the Born level in the electroweak interactions. 
We take the m, d, s, c and h quarks to be exactly massless. We include contributions from 
virtual top quarks that behave like log"' {mf / M^) , dropping terms that behave like {Mz/mf)" 
as rrit oo. 

Given this theoretical framework, the theoretical expression for R{M'^) has the form 



Here Rq is the value of R in the parton model, without perturbative QCD corrections. The 
QCD corrections are contained in 7^(M|) = as(M|)/7r + 0{al), which is often denoted 
(^QCD- We study 7?.(M|) and try to estimate the theoretical uncertainty in 7?.(M|) caused 
by evaluating it in perturbation theory truncated at order a^. For this purpose, we use a 
nominal value asiM^) = 0.120 of the MS strong coupling evaluated at Mz- If an experimen- 
tal value for 7^(M|) were used to extract as{M'^), then the fractional theoretical uncertainty 
in 7^(M|) would translate into a fractional uncertainty of the same size for a^(M|). 

When we present numerical results, we choose Mz = 91.188 GeV and sin^ 9w = 0.2319. 
We take the top quark pole mass to be 170 GeV, as estimated in Ref. ||l| from the CDF and 
DO resuhs §. 

The scope of this paper is limited, and in fact we do not attempt to evaluate 7^(M|) at 
the level of precision that we are discussing. Such an evaluation involves careful consideration 
of a large number of small effects. Among these are electroweak effects beyond the Born 
level 0, effects of non-zero masses for the light quarks and (M^/mf)"' contributions from 
virtual loops containing the top quark [^j. We review the status of some of these issues in 
the appendix to this paper. 
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II. THE RUNNING COUPLING AND TOP MASS 



In this paper we denote by as{s) the running couphng in a renormahzation scheme that 
may or may not be the MS scheme We denote by Cis{s) the running couphng as defined 
by the MS scheme with five fiavors of hght quarks. 

The dependence of as{se^) on t is given by the renormahzation group equation 



T+ i 1 tj = JT^ /5(«s(se*)) 

at \as\se') ) a;;(se*) 




= /3o + /3i^^+/32M^ +■•■• (2) 

TT y TT 

We use this equation to derive an approximation for 7r/Q;s(se*). We find 

+ + ^log(l + a;) 







, «.(g) / /9o/32-/3? X ^ log(l + x) 



+ ■■■, (3) 



where 



x = /3ot^. (4) 

TT 

Further terms in this series involve higher powers of as(s)/7r times functions of x that 
are proportional to x for small x. We do not include any more terms because the next 
term involves the coefficient /^s, which is unknown. If we wanted to recover the ordinary 
perturbative expansion of 7r/as(se*) up to order ^^(s)^, we would note that x is proportional 
to as and expand in powers of x, then omit terms beyond x^ or xas- Eq. is better than 
the purely perturbative expansion because it is a valid expansion in powers of asi^s) when x 
is fixed at some finite value. Thus it is useful when (i^t is as large as 7r/as(s). 

We shall sometimes want to examine the dependence of the results of calculations on the 
renormalization scheme used in the calculation ((7/. Ref. |]^). For this purpose, we define 
an as in a renormalization scheme that may not be the MS scheme by 

as{s) = as{s) + C2as{sf + c^asisf H . (5) 

Then one can use as{s) as the expansion parameter of the theory. Since the perturbative 
formulas used are inevitably truncated at some order of perturbation theory, the results 
depend on the coefficients q that specify the scheme. We will want to find out how much the 
results depend on the q. There are two purposes to this. First, the choice of renormalization 
scheme represents an ambiguity of the theory, and we want to have an estimate of the 
numerical importance of this ambiguity. Second, there are uncalculated higher order terms 
that are, by necessity, omitted from the calculation. Parts of these terms serve to cancel the 
dependence of the results on the q. Thus the observed size of the dependence of the result 
on the Ci serves as a rough indicator of the size of the uncalculated higher order terms. 
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The coefficient C2 can be simply absorbed into a change of the scale of the running 
coupling: 

a,{s) = as{se^^) + 4 a^(se^*) + • ■ ■ . (6) 

That is, using Eq. @ on the right hand side of Eq. reproduces Eq. (^). The term 
in Eq. (||) proportional to results in a change of the coefficient /?2 in the P function 
that describes the running of a^. (Recall that Po and Pi are scheme independent.) Let us 
parameterize this change as 

P2 = P2 + SP2, (7) 

where P2 is the third coefficient of the P function in the MS scheme and other MS-type 
schemes. Then the relation between as and as can be written as 

We shall use 6t and 6P2 to parameterize the choice of scheme. 

By combining Eq. with Eq. (]^), we see that as{s) can be expanded in terms of as{s) 
by using 

^ = ^(l + 5^) + |log(l + 5a;) 
as{s) as{s) Po 

^ Ms) ( P0P2-I3! Sx ^ Pi log(l + 5x) _ 5P2 1 \ 

TT V P^ 1 + Sx Pi l+5x Po l + 5x) 

+ ■■■. (9) 



Here, 



5^ = /3o^t^. (10) 

TT 



In the framework of this paper (except for the appendix), light quark masses do not 
appear in R{s) because they are set to zero. However, the top quark mass does appear, 
starting at order a^. Thus it is necessary to state carefully how we define m^. We let mt{s) 
be the running top quark mass within the MS scheme. At the level of perturbation theory 
at which we work, we need the one loop evolution of mt(s), which we write as 



« m?(M|) exp ( 2.„ r '4 ^) « ™|(M|) ) ' (11) 



270 /A) 



with 



70 = -1. (12) 



(See, for instance, Ref. |^). One can, of course, use a different scheme and define a running 
mass 
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m2(s) = mf{s) [1 + Cias{s) + ■■■]. 



(13) 



We do so, absorbing the first coefficient Ci into a change of scale by an amount 61^. Thus 
we define 

(— / \ \ 27o//3o 

The parameter St^ can be chosen independently from the scaling parameter 5t in the defi- 
nition (H) of the coupling. 

The dependence of TZ{s) on the top quark mass is quite small, so the dependence of 
7?.(s) on 6tm is also small. In fact, we find that 7^(M|) varies by only 0.3 parts per mill 
for —4 < 6tm < 4. In order to limit the parameter space to be explored in our numerical 
examples, we therefore set 

Stm = 0. (15) 

Thus the running top mass at s = M| in our examples is simply the MS running top mass 
mt(M|). We take mf (M|) = 170.2 GeV, which corresponds to a pole mass of fht = 170 GeV 
after use of 

TJit [Mz) ^ exp ( -- — - — ) ( ^JJJT^ ) ■ (16) 



(See, for instance, Ref. p). The value 170 GeV is estimated in Ref. ||ll| from the CDF and 
DO resuhs 0. 



III. PERTURBATIVE EXPANSIONS 

With the theoretical framework defined in Sec. |, the theoretical expression for R{s) has 
the form 

R{s) = R^{l + n{s)}. (17) 

Here Rq is the value of R{s) in the parton model, without perturbative QCD corrections. 
The QCD corrections are contained in TZ{s), 




7^(s) = 7^l ^ + 7^2(^^) +7^3(^) +■■■. (is) 

The value of 7^(s) calculated in finite order perturbation theory depends on the parameters 
6t, 5(32 and 5tm that define the renormalization scheme. We keep these parameters arbitrary 
in this analysis in order to be able to test the sensitivity of the calculated value of 7?.(s) to 
their choice. As already noted, the dependence of 7^(s) on 5tm is negligible. 

The t dependence of a<i(se*) is given by the renormalization group equation (|^). The 
coefficients of the /? function that appears in this equation are 



5 



/5o = (33-2iV^)/12, 
(3^ = (306 - 38A^/)/48, 

p2 = (77139 - 15099 Nf + 325 ArJ)/3456 + Sp2 , 

where Nf = 5 is the number of hght quark flavors used throughout this paper. 
The coefficients 7^1,7^2, "^3 are [jl^'O 



(19) 



7^l 



^-llC(3)-iV, 
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2C(3) 



1 
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, Ml 



87029 1103 C(3) 275 C(5) /^qV 



288 
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7847 262 C(3) 25C(5) 
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5C(3) 



-18.65440 
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31 



151 19C(3) 



162 



27 



log I 



+/?2(<5t)2 + [A + 2/?o7e2,o] 5t-^-^. 

Po 



. Ml 



23 , 2, 
+ ^log^l 



(20) 



+ 270 5tr 



Here (f j, a^) with i = u,d, s, c, 6 is the (vector, axial vector) coupling of the quark of ffavor 
i to the Z boson, as specifled in the appendix. We use 7^2,o to denote 7^.2 with 6t = 0. We 
recall that the mass anomalous dimension is 70 = —1 and that the parameters 6t, 5(32 and 
5tm give the scheme dependence, as described in Sec. |I[ 
The numerical values are 



/?o ~ 1.92, 
(3i ^ 2.42, 
(32 ~ 2.83 + 6(32, 



(21) 



and (with St^ = and mt{s) = mt{Ml)) 



7^l = l, 

7^2 ^ 0.76 + 1.92 

7^3 ^ -15.73 + 5.35 St + 3.67 (StY - 0.52 5/52. 



(22) 



In this paper, we will deflne various approximations TZa to TZ. The flrst of these is the 
simple third order perturbative approximation: 



7^^(.;Pert) = E7^,(^^y 



(23) 



6 



As discussed in Sec. ||, the renormalization scheme ambiguity can provide an estimate, or 
at least a lower bound, on the theoretical uncertainty produced by truncating perturbation 
theory at order a^. To investigate this ambiguity, we show in Fig. |1| a contour plot of 
7?.^(M|; Pert) as a function of St and SP2, with 6tm = 0. The range shown for the scale 
parameter, —4 < 6t < A, corresponds to scales fi = [se^]^^"^ in the range O.IA Mz ~ M^e~^ < 

I. The range shown for 6/32 corresponds to schemes with 



fi < Mze^ ^ 7AMz in Eq. 



6.8. 




j_i_uL 
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FIG. 1. Contour plot of the simple third order perturbative approximant 7^yi(M^; Pert), 
Eq. (p3|), versus the scheme fixing parameters 6t and 6(32, with 6tm = 0. 



We learn from Fig. |I| that TZ aIM"^; Pert) is not very sensitive to 5/?2- Accordingly, 
we set 6(32 = and plot 7^A(M|; Pert) versus 6t in Fig. |. We note that 7^A(M|;Pert) 
varies by about 2.6% between its local maximum and its local minimum. We conclude 
that TZ{M^) probably lies within this 2.6% range. Thus we ascribe a theoretical error of 
±2.6%/2 = ±1.3% to the value of 7^(M|). In the remainder of this paper, we attempt to 
reduce this error by using more sophisticated methods than simply taking the first three 
perturbative terms in 7?.(M|). 

The perturbative series for 7^(M|) provides our starting point. We see that the series 
for (3{as) is nicely behaved, but that the series for 7l{s) not as well behaved, with a large 
value for TZ^ at St = 5(32 = 5tm = 0. In fact, this large value can be attributed to the term 
-/?^7rV3 ^ -12.1. 
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FIG. 2. Plot of the usual third order perturbative approximant 7^ /^(M^; Pert) versus the 
scheme fixing parameter 6t with 6P2 = and 6tra = 0. 



IV. 7r2 TERMS 



The offending vr^ term in TZ^ arises, at a rather mechanical calculational level, because 
factors of log(— s ± ie) = log(s) ± in occur in the calculation, leading to powers of vr in the 
result. In order to see what happens at higher orders of perturbation theory, we write R{s) 
as a discontinuity: 

R{s) = — |n(-s + ie, /x^) - n(-s - ie, . (24) 

Here C is a normalization constant and Il{Q^, fi'^) is the standard Z boson self-energy func- 
tion including the QCD contribution. It is proportional to the Fourier transform of the time 
ordered product of two weak current operators. The current operators carry momentum q^. 
We define = —q^q^, so that > if the momentum q^^ is spacelike. The function 11 
depends on the renormalization scale n^. However the function 

^(Q') = -Q'^^(Q^/.2) (25) 

is a renormalization group invariant. The derivative here avoids the overall renormalization 
in H. For this reason, it is standard practice to work with D{Q'^). 
We may write the perturbative expansion of D{Q'^) in the form 

D{Q^)=D,[l+V{Q^)], (26) 
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where Dq is the value of D in the parton model and where 



asm 



+ V2 



asm 



asm 



+ ■■•. 



(27) 



■1 



The first three coefficients T>n are the same as the corresponding in Eq. ( pOD if one 
substitutes for s, except that lacks the term — /^qVt^/S. The numerical values (with 
5tm = and mt{s) = mi(M|)) are 



If we stay near 6t = 0, this series appears to be quite nicely behaved. We believe on the 
basis of general arguments (to be discussed in the next section) that the coefficients Vn will 
eventually grow for large n. However, that growth is not apparent in the first three terms. 

The function D{Q'^) is calculated using Euclidean quantum field theory, in which only 
very weak infrared singularities occur near the contour of the internal momentum inte- 
grations. On the other hand, a direct calculation of R{s) involves Minkowski momentum 
integrations over regions in which various internal particles can go on shell. Only some deli- 
cate cancellations prevent R{s) from being infinite. Surely D{Q^) should be better behaved 
than R{s). This observation leads to the following 

Hypothesis 1. The perturbative expansion of D{Q'^) remains well behaved beyond 
the three terms that are known, subject only to the eventual growth of the Vn 
dictated by the standard renormalon and instanton ideas. 

We adopt this hypothesis here, although it is criticised in Ref. [|1^ on the grounds that there 
could be other sources of large perturbative coefficients in T'(s). 

We are interested in the observable function R{s). If we accept this Hypothesis 1, then 
instead of calculating R{s) directly, we should relate it to the nicely behaved function D{Q'^). 
From Eqs. (|2^) and (^) we obtain 



In the following section, we will deal with the expected large order behavior of the T>n by 
following the standard practice of using the Borel transform T) of T): 



Pi = l, 

V2 ^ 0.76 + 1.92 5t, 

V3 ^ -3.65 + 5.35 5t + 3.67 {5tf - 0.52 5p2- 



(28) 




(29) 




(30) 



If we write the perturbative expansion of T^lz) as 



V{z) = Vo + Vi Z + V2 + ■■■, 



(31) 



then 
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(32) 



Because of the l/n\ factor, the perturbative expansion of ID in powers of z is much nicer 
than that of T) in powers of O-sji:. In fact, one expects T'i^z) to be analytic near z = 0. As 
discussed, for example, in Ref. |]T3, there are singularities expected in the complex z plane. 



including some on the integration contour along the positive z-axis. In addition T>{^z) is not 
expected to be well behaved as z ^ oo. Thus the meaning of the integration in Eq. ( pUD is 
ambiguous. In this section, we simply leave it as ambiguous. 
We can relate 7^(s) to V by inserting Eq. ^ into Eq. (|2|): 



oo 



where 



with 



n{s)= I dz expi ^] F{asis),z)V{z), (33) 

Jo \ aJs] 



F(a„ z) = ^ r de exp(-zG'(a„ 9)) , (34) 



Giasis), 6) = — (35) 

Eq. (|33D is the basis for the analysis in this paper. We note that the factor 7c/as{s) in 
the exponent in Eq. (^) is big, about 30 for ~ 0.12. Therefore the integral over z is 
dominated by small z, PqZ < Poag/n ~ 0.06. Thus we will be primarily concerned with the 
expansion of 'Diz) in powers of z. 

Before addressing T>{z), however, we need a good approximation for F{as,z). Since 
small z is important, we are particularly interested in the small z region. However, it 
is rather easy to find an approximation for F{as, z) that is good for a wide range of z, 
based on the smallness of its argument a.^- We use the solution (|^) of the renormalization 
group equation (|]) for vr/a^ in order to derive an approximation for G{as,0). We find 
G{as,0) ^ GA{as,0) where 



GA{as,e) = t(3o0 + ^iog{l + y) 




where 



y = iPoe-. (37) 

TT 



Further terms in this series involve higher powers of cts/vr times functions of y that vanish for 
?/ — s> 0. The ordinary perturbative expansion of G{s) results from expanding in powers of y, 
which is proportional to a^, then omitting terms beyond y^ or yas- However, as{s)/Ti ~ 1/30 
while \y\ 1/5 at ^ = vr. Thus as{s)/TT is a much better expansion parameter than y. Since 
we don't have to expand in y, we don't. 
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We now have an approximation Gyi(as, 6) for G{as, 0). Our corresponding approximation 

FA{as,z) for F{as,z) is 



(38) 



FA{as,z) = — d6 exp{-zGAias,9)) , 



with the integral computed to sufficient accuracy by numerical methods. In Fig. |] we show a 
graph of FA^as, z) versus (3qZ superimposed on a graph of exp(— vr^/a^), all with = 0.12. 



. 8 
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-0.2 




exp(-7r2^/Q:s) 
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1 . 

(3qZ 



1 . 5 



FIG. 3. Graph of Fa^oisiZ) and exp(— vrz/as) versus j3oz with = 0.12 

How good is our approximation ^^(as, z)l The first omitted term in G{as-, t 
MS scheme, 



is, in the 



where 



h{x) 



AG{as) = ( ^ ) h{x) 



Pf log'(l + x) , pip2 log(l + X 



(39) 



+ 



1 + x)^ ' (3i {1 + xy 

x^ P1P2 X , Ps x{x + 2) 



+ 



2Pl{l + xy Pi (l + x) 2Po{l + x) 



(40) 



This term contains a factor (a^/Tr)^ 10"'^ for ^ 0.1. This factor multiplies h{x), which 
cannot be evaluated because it contains the unknown coefficient P3. However, we can see 
from the structure of h{x) that it is not large unless /?3 is large. We can get a quantitative 
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idea of the effect of AG by choosing some plausible values for P3 and then calculating 7l{s) 
with AG included. We find that, taking (3^ in the range —10 < /^s < 10, the fractional 
change 7?.(M|) induced by including AG{as) is no larger than 2 x 10~^. Since this error is 
small compared to our target error of a few per mill, we can safely neglect it. 

We thus obtain an approximation for TZ that uses third order perturbation theory but 
sums certain "vr^" effects to all orders: 

roc I -j^^ \ _ 

nA{s-Ti^)= dz —] FA{as{s),z)VA{z), (41) 

Jq \ as{s) J 

where Fa is given in Eq. (|38|) and Va{z) is simply T>{z), Eq. (|31|), expanded to second order 
in z. 

This treatment of tt^ terms is similar in spirit to that of Le Diberder and Pich , who 

) in perturbation theory, use Eq. (|]) for ^^(se*^), and perform the 
9 integral exactly. We simply embed this approach into the Borel transform. 

In Fig. I we plot Te^lMf; TT^) versus the scheme parameter 5t with the other scheme 
parameters set to 5(32 = and 5tm = 0. We overlay the plot of 7^^(M|; Pert) from Fig. |^. 
We note that TZAiM^; TT^) varies by about 0.32% between its local maximum and its local 
minimum. This is a much smaller variation than that of 7^^(M|; Pert). A very optimistic 
view would be that 7^(M|) probably lies within this 0.32% range, so that one would ascribe 
a theoretical error of ±0.32%/2 = ±0.16% to the value of 7^(M|). However, this error 
estimate is smaller than other error estimates that we will develop later. We therefore 
regard the flatness of the curve for 7^^(M|; vr^) as being partially the result of an accidental 
cancellation, and refrain from taking 0.16% as a reasonable error estimate. 

We close this section by emphasizing the observation that the straightforward pertur- 
bative expansion of 7^(s) is, in part, an expansion in powers of (3o0 [«s('S)/vr], with 6 ^ tt, 
instead of an expansion in powers of [as(s)/7r] ^ 1/30. One can attribute the appearance 
of "vr^" terms in TZ{s) to this phenomenon. This observation helps to make Hypothesis 1 
plausible. Unfortunately, this argument is only suggestive, since one can not be sure that 
there are not "bad expansion parameters" lurking somewhere in the calculation of V{Q^). In 
the next section, we turn to the behavior of the perturbative coefficients in 'Dlz), assuming 
that the evidence for a bad expansion parameter is not, in fact, lurking just beyond the last 
calculated coefficient. 

V. TRUNCATION OF THE INTEGRAL 

If we do not expand 'Dlz) in powers of z, then at this point we have an approximation 
for 7?. of the form 

/■oo / jYZ \ 

7^(s) ^ / dz exp — FAiasis),z) V{z), (42) 

with Fa given in Eq. (PH]). Since as{s) is small, the dominant integration region is z <^ 1. 
Indeed, taking ~ 0.12, we have exp^—nz/as) < 10~^ for Pqz > 0.51. Thus it is useful to 
write TZ in the form 



expand V{se''^) in Eq. (29 
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FIG. 4. Plot of TT^-summed approximant TIa{M'^;tt'^), Eq. (pl|), (dashed line), versus the 
scheme fixing parameter 6t with 6P2 = and 6tm = 0. We also show 7?.yi(M^; Pert) from Fig. ^ 
(full line). 



7^(s) ^ / dz exp — Fa(«.(s), z) V{z) + 7^/^(s) (43) 

Jo \ as[s) ) 

with /?o-2^max ~ 0.5. The fundamental question of how the "sum of perturbation theory" is 
precisely defined relates to the definition of TZr. In turn, this question is related to how the 
renormalon and instanton singularities are treated and to the question of the convergence 
of the integral at large z. However, our purpose here is at once more modest and more 
practical. We adopt 

Hypothesis 2. It is safe to ignore the large z part of the Borel integral when 
calculating 'R.{s) for s ~ M^, even though this part of the integral is ill defined. 

We thus neglect TZr and concentrate on the integral up to Zmax in Eq. (|i3|). The advantage 
is that we can use approximations for ^^{z) that have singularities on the positive z axis 
outside of the region of integration. 

We can test the sensitivity of the computed value of TZ to z^^^y, by replacing V{z) by its 
second order expansion in powers of z. Then the ratio of the two terms in Eq. (^) with 
^max = 0.5 (for s = Ml, 5t = 6t^ = 6(3s = 0) is 7^/^(s)/7^A(s; tt^) ^ 6 x 10"^. 
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VI. ACCOUNTING FOR RENORMALONS 



We now turn to the perturbative expansion 



(44) 



The coefficients Vn can be expressed as an integral, 



/ dz z-^'-^Viz). 

2-iTi Jc 



(45) 



The contour C encloses the point z = but excludes any singularities of T^lz). Thus the 
behavior of the P„ at large order n is controlled by the part of the contour that lies nearest 
to 2 = 0, which in turn is controlled by the singularities of 'Diz) that are nearest to z = 0. 
A singularity of the form (, 



n 



A-l 



- zq)~ makes a contribution to P„ that is proportional to 
Thus the most important determinant of the singularity's contribution to the Vn 
at large n is its location, zq. A small zq produces large coefficients. The next most important 
determinant is the strength of the singularity, A. A large positive value of A produces large 
coefficients. 

The nearest singularities are thought to be the first two ultraviolet renormalon singular- 
ities at Poz = — 1 and —2 and the ffist infrared renormalon singularity at (3oz = +2 ||13|,|15 



In this section, we use the available information on these singularities to obtain a pertur- 
bative expansion that has better convergence properties. Of course, "better convergence 
properties" refers to the perturbative coefficients for large n. It is problematical whether 
convergence improvement helps already after only three terms of the series. 

The ffist ultraviolet renormalon singularity is at Pqz = —1. This is the singularity that 
is closest to the origin (at least so far as anyone knows). It thus controls the large order 
behavior of the perturbative series. Unfortunately, the theory of the ultraviolet renormalon 
singularities is not as simple or as well developed as that for the infrared renormalon sin- 
gularities. (See, however, Ref. ||16[). For instance, the strength of the singularity is not 
known. 

The ffist infrared renormalon singularity is at jSoz = +2. There are other singularities 
farther away from the origin along the positive real z-axis, but we need not be concerned 
with them: since they lie farther from z = 0, their contribution to the large order behavior 
of the perturbative coefficients is weaker than that of the ffist singularity. It is significant 
that there is no infrared renormalon singularity at /SqZ = +1. The ffist infrared renormalon 
singularity has a power behavior. 



Viz) 



Poz 
2 



1 -1 



-2/3i//3g 



(46) 



where c is a constant [|1^,|T^. Numerically, the exponent is —1 — 2(3i/ (3q 



-2.3. 



We can make use of this information. Consider the function 



C{z) = V{z) 



2 



l+2/3i//32 



(47) 
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The factor multiplying T>{z) cancels its divergence as Pqz 2. The function C{z) is still 
singular at (3qZ = 2, since if we multiply a term in V^z) that is analytic at PqZ = 2 by the 
nonanalytic factor, we create a nonanalytic term. However, singularity is much weaker than 
it was, behaving like 



Ciz) 



2 



l+2/3i//3g 

(48) 



Thus the perturbative expansion of C{z) would be better behaved than that of 'Diz) at 
large order if it were not for the fact that the leading ultraviolet renormalon singularity at 
PqZ = — 1 dominates the large order behavior. 

We can, however, improve the large order behavior arising from the leading ultraviolet 
renormalon by merely moving it out of the way by means of a good choice of variable. 
Following Mueller [T^ we define a new variable ( by 



= T\ /-IA\2 ' /^oC = 4 , -■ (49) 

(1 - M/^r VI + Pqz + 1 

This transformation maps the origin of the 2;-plane onto the origin of the C-plane. We have 
chosen the normalization of C such that 

Cr^z + 0{z^). (50) 

near z = Q. The map treats specially the interval jd^z < — 1 on the negative z-axis that 
contains the ultraviolet renormalon singularities. The whole complex 2;-plane except for this 
interval is mapped to the interior of the disk |/?oCI < 4 in the ^-plane. The singularity- free 
interval — 1 < fS^z < in the negative z-axis is mapped onto the interval —4 < /5oC < of 
the negative (^-axis while the interval < Pqz < oo on the positive z-axis, which contains 
the infrared renormalon and instanton singularities, is mapped into the interval < /5oC < 4 
of the positive ^-axis. 

We consider the function 

m = c{z{0). (51) 

The singularity of B{Q that is nearest to the origin of the C-pla^ne is the first infrared 
renormalon singularity, which is at 

/3oC = 4^^^^ ^ 1.1. (52) 

Thus moving the ultraviolet renormalon singularity away has had a price. We have moved the 
infrared renormalon singularity closer to the origin. However, we have previously softened 
the infrared renormalon singularity, so the price is not too great. The net effect should be 
an improvement. 



The effect of singularity mapping has been investigated recently by Altarelli et al. [|T2 
However, these authors did not also soften the infrared renormalon singularity. They found 
that there was no gain in this method. 
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In order to use the singularity softening and mapping, we use the first three terms in the 
perturbative expansion of V^z), 



V{z) = 1 + 

+ 



0.40 + 1.00 (5t 



-0.50 + 0.73 (5t + 0.50 {5tf 



(53) 



to calculate the first three terms in the perturbative expansion of B{Q. The result is 



-0.76 + 1.00 5t 



+ 



-0.96 + 0.07 5t + 0.50 {5tf 



(54) 



(Here we have displayed the coefficients numerically, with the choices 6(32 = 6tm = and 
mt{s)=mt{Ml).) 

This perturbative series for B{() is supposed to be better behaved at large orders than 
was the perturbative series for T>{z). The expected improvement is not, however, visible in 
the first three terms. In fact, we started with a series that was quite well behaved, and we 
have applied a rather mild improvement program. As long as the infrared and untraviolet 
renormalon singularities are as described in this section, this program may be expected to 
make the perturbative coefficients smaller at high order, but one cannot expect too much 
to happen at order two. 

An example of this procedure applied to a simple model may be useful as an illustration 
of what happens at high order. Suppose that 



V{z) 



[l + M [l-(3oz/2f 



(55) 



with p 



1 + 2/3i//5o. Then the perturbative expansion of T>{z) is 

V{z) = 1 + 1.68 /3o^ + 0.44 (/5o^)' + 1.21 {(3ozf 

-0.06 {(3oz)^ + 0.81 {(3ozf - 0.35 {(3ozf + 



(56) 



Applying the renormalon improvement procedure gives the function B{() with a perturbative 
expansion 



BiC) = 1 + 0.52 (3oC - 0.87 {M? + 0.30 (M? 

-0.02 (/?oC)' + 0.06 iPoC? - 0.01 iPoCf + 



(57) 



The series for B is clearly better behaved at high orders than the series for P . One might 
claim to see an improvement beginning with the fourth term, which corresponds to the first 
uncalculated term in the case of the real V and B functions. However, at this quite low 
order of expansion, the improvement is marginal. 

The procedure for singularity softening and mapping may be summarized as follows. We 
calculate the first terms in the expansion of B{() according to Eqs. (0) and (pi]), where 
for us A^ = 3. Then we instead of using 
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n=0 



for V{z) in Eq. (^31), we use 



/3oz 



n=0 



This gives an approximation for TZ{s) that we may call TZAis; vr^, R'lons): 



7^yi(s; vr , R'lons) 



(i^; exp 
\ as[s) 



FA(aJs), z) 



[l-Poz/2] 



l+2/3i//3^ 



(58) 



(59) 



(60) 



The replacement of T>a{z) by Va'{z) does not modify the integrand much. In Fig. ^, we 
show the ratio Va'{z)/T>a{z) as a function of Pqz. We see that this ratio is nearly 1.0 in the 
important integration region (3qZ < 0.2. 



'6 



0.9 



0.967 




0.1 0.2 0.3 0.4 

FIG. 5. Modification of the Borel integrand to account for renormalons. We plot 
'Da'{z)/'Da{z), Eqs. ( ^8|) and (^), versus Pqz. We set s = M| and choose the scheme fixing 
parameters 6t = 6^2 = ^tm = 0- 



VII. RESULTS 

We have developed an approximation to TZ{s) that takes vr^ contributions into account 
and uses information about the leading renormalon singularities to try to improve the con- 
vergence of the perturbative expansion for 'Di^z). In Fig. ^ we plot this approximation, 
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versus the scheme parameter St with the other scheme parameters set to 
Sj32 = Stm = 0. We overlay the plots of the pure perturbative function, 7?.^(M|; Pert), and 
the approximation that simply takes vr^ contributions into account, 7^A(^I;7r^). We note 
that 7^A(s;7^2,R'lons) varies by about 0.8% between its local maximum and its local min- 
imum. This suggests that 7?.(M|) probably lies within this 0.8% range, so that one would 
ascribe a theoretical error of ±0.8%/2 = ±0.4% to the value of 7?.(M|). 



. 039 



. 03 



. 037 




FIG. 6. Plot of approximant 7^a(s; ^r^, R'lons), Eq. (|60|), versus the scheme fixing parameter 
6t with 6P2 = and Stm = (heavy line). We also show 7^^(M|; Pert) from Fig. |2| (light line) and 
71^(^1; 71-2) from Fig. | (dashed line). 

We can take another approach to error estimation. We note that the first three co- 
efficients of (/5oC)" in Eq- (HD are all of order 1. That the coefficients do not appear to 
be growing or shrinking with n is normal since the series is expected to have a radius 
of convergence of about 1 in the variable Po(. We thus expect that the uncalculated co- 
efficient of {PoCY will also be of order 1. If we add a term 1 x {PoCY to the series in 
Eq. (0), 7?.^(M|) changes by an amount STZa that can serve as an error estimate. We find 
SnA/nA{Ml- TT^, R'lons) ^ 0.2%. 

We thus have three error estimates. From the St dependence of 7^yi(M|; vr^) we estimated 
a 0.16% error. From consideration of the likely size of the next term in B{C,) we estimated a 
0.2% error. From the St dependence of 7^^ (M|; vr^, R'lons) we estimated a 0.4% error. We 
take the largest of these values, 0.4%, as a reasonable estimate of the theoretical error (in 
the spirit of a "1 a" error). 

For the central value, we take the value of 7^a('5; tt^, R'lons) at St = 0, which is almost 
exactly also the value of TZa{s] ti"^) at St = 0. This value is 



7^A(M|;7^^ R'lons) 



st=o 



;i - 0.006) X 7^A(M|;Pert) 



<5t=0- 



(61) 
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That is, our best estimate for TZ is renormalized down by 0.6% compared to the standard 
MS value with a scale choice /i = Mz- 

One often uses a measurement of 7?.(M|) to extract a value of as{M'^). Recall that, 
to a good approximation, 7^(M|) oc as{M'^). Thus the value of ag^M"^) extracted from 
data using the "standard" MS expression for 71 (with a scale choice /i = Mz) would be 
renormalized up by 0.6% if one uses the "improved" version of TZ presented here: 



[as(M^)]improved ~ 1-006 [as(M^)]standard- (62) 

The fractional error to be ascribed to as{M^) from uncertainties in the QCD perturbation 
theory is just the fractional error in 7?.a(M|; tt^, R'lons) estimated above as 0.4%. This is 
one third of the 1.3% error that we would ascribe to as(M|) extracted using the standard 
perturbative approximant 7?.^(M|; Pert). The shift in Eq. (|6^) is about the same size as the 
estimated theoretical error, so it is marginally significant. 

We note that the experimental error for the extraction of as by this method is about 
5% [|T7|, much larger than the QCD theoretical error that we estimate above. There are 
also sources of theoretical error not associated with QCD. According to the estimates of 
Hebbeker, Martinez, Passarino and Quast [jl8[, the most important of these are a ±2% 



uncertainly from electroweak corrections and a ±2% uncertainly from not knowing the Higgs 
boson mass. 
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APPENDIX: PRESENT STATUS OF PERTURBATIVE QCD EVALUATION OF 

Z DECAY RATES 



The decay rate of the Z boson into quark antiquark pair can be written in the following 
form: 



Z— »hadrons 



8V27r 



+Pfa} 



1 _ 4X^)3/2 + J^^^ Xf, Xt) + Cd(«, X 



qedI«,«s,-^/J 



(Al) 



Here there is a sum over light quark flavors f = u,d, s, c, b. We define Xf = 'mf{Mz)'^ /M"^ 
and Xt = mt{MzY /M"^. (We use the MS definition of masses.) 
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The vector and axial couplings of quark / to the Z boson are f / = 

and CLf = 21 J . The electroweak self-energy and vertex corrections are absorbed in the 
factors pf and kf. The current status of the electroweak contributions has been discussed 
in detail in Ref. Q]. The small QED corrections in vector and axial channels have the form 

Co = 7 4 ? [1 + 12^/ + + + 0{aas\ (A2) 



"qed 



- e) - [1 - QXf - 12Xf logXf + 0{Xj)] + 0{a^) + 0(aa,). 

4 TT ■' 



(A3) 



The corrections of order and aa^ are discussed in Ref. . 



It is convenient to decompose the QCD contributions into singlet and non-singlet parts 
and further into vector {V) and axial vector {A) contributions. The nonsinglet parts are 
represented by the terms 6^,;,^ and 5qcd) correspond to cut Feynman graphs in which 
a single quark loop of flavor / connects the two electroweak current operators. The singlet 
contributions correspond to graphs with the electroweak currents in separate quark loops 
mediated by gluonic states. In the singlet contributions one does not have a single sum over 
a flavor /. These contributions are represented by the terms and C^. 

The nonsinglet QCD contribution in the vector channel to order can be written in 
the form 



+ 



a. 



TC 



1.40923 + 104.833 + Y,F^^\Xv) + G^'^\Xt 
-12.76706 + 547.879 + Y,F^^\X^) + C^^H^t) 



(A4) 



In this formula, as denotes the running MS coupling in five flavor theory evaluated at Mz- 
The transformation relation for different number of flavors and different scales, as well as 
the relation between the MS running mass and the pole mass can be found in Ref. ||20[ . 



The order and a^ terms have been evaluated in the limit of vanishing light quark 



masses and infinitely large top mass in Refs. |jTO|,[n]|. These contributions, (a^/Tr) + 
1.40923 (as/vr)^ — 12.76706 (a^/Tr)^, are the {vector, nonsinglet} part of the perturbative 
series analyzed in the main body of this paper. 

The terms proportional to Xf represent the leading corrections to the approximation 
Xf = 0, as given in Ref. 0]. 

The function F^^\X^) arises from three-loop diagrams containing an internal quark loop 
with a quark of flavor v = u,d, s, c, b propagating in it (while the quark of flavor / couples to 
the weak currents). This function represents the corrections to the approximation X^, = 0. 
These contributions are already small, so it suffices to approximate Xf by in F^'^^X,^). In 
fact, numerically 0, 



(A5) 



F(2)(XJ ^ X^ X I -0.474894 - logX^ + [-0.5324 + 0.0185 log X^]| 
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is so small that the whole function could be neglected. 

The function G^'^\Xt) represents the contribution of virtual top quark loops inside three- 
loop cut Feynman diagrams. These contributions are small since the top quark is nearly 
decoupled from the theory. Thus it suffices to approximate Xf by in G^'^\Xt). Numerically, 
one finds 



r 44 
1675 



+ 



135 



\ogXt - V^r^ [0.001226 + 0.001129 log 



(A6) 



The first two terms in the right hand side of eg . ( |A6|) have also been obtained using the large 



mass expansion method |^ . 

At order a^, there can be two internal quark loops. However, it suffices to consider only 
one loop with a nonzero light quark mass at a time, or one top quark loop with all light 
quark masses set to zero. Then we can define functions F^'^^X^) and G^^\Xt) analogously 
to and G^^\Xt). For F^^\X^) the following small mass expansion is obtained in 

Ref. i: 

^ -6.12623Xj. (A7) 

For G^^\Xt) the following large mass expansion has been obtained in Ref [22|: 

G^^\Xt) ^ Xf^ X [-0.1737 - 0.21241ogXt - 0.0372 log^ X^]. (A8) 

The nonsinglet contribution in the axial channel is the same as the one in the vector 
channel except that the contributions proportional to Xf [Q,^ are different: 



rA.ns 
"qcd 



a 



TT 
+ 



^ [1-22X^] 

1.40923 - 85.7136 Xf + Y. + C^^H^t 

V 

-12.76706 + (unknown) Xf + J2 F'-^HX^) + G^^\Xt 

V 

We now turn to the singlet contributions, which start at order a^: 





)1 















(A9) 



a 



TT 



,V/A 



a. 



TT 



+ 



At order a?, there is no vector contribution, 



>~'2 



0, 



(AlO) 



(All) 



while the axial contributions from u and d quarks and from c and s quarks vanish in the 
limit of vanishing quark masses. This is because in the Standard Model the quarks in a weak 
doublet couple with the opposite sign to the weak axial current. However, the contribution 
from the t,6 doublet is significant because of the large mass splitting 113: 
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37 7 

= ~ 12 ~ + 81^*"' + O.OlSXr' 

Xft /80 



+X, (18 + 6 logX,) - ^ + _ logX, ) . 



(A12) 



Here the corrections proportional to X;, have been calculated in Ref. 



At order a^, both channels contribute. The vector contribution in the limit of massless 



light quarks is [11 



£y = -0.41318 (^w/j + (0.02703 + 0.00364 X-2 + 0(X-=^)) (A13) 

The sums here run over light quark flavors f = u,d, s, c, b. The terms proportional Xf^,X^'^ 
were computed in Ref. |2^ and turn out to be negligible. 

In the axial channel, the order singlet contribution in the large top mass expansion 



reads [25,22 



^3 



-15.98773 - ^ log + ^ log^ X^. 



(AM) 



Corrections for a nonzero b quark mass are not yet known. However, at the level of precision 
of this paper, they are not expected to be significant. 
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